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We study the duality of reconstruction systems, which are g-frames in a finite di- 
mensional setting. These systems allow redundant linear encoding-decoding schemes 
implemented by the so-called dual reconstruction systems. We are particularly inter- 
ested in the projective reconstruction systems that are the analogue of fusion frames 
in this context. Thus, we focus on dual systems of a fixed projective system that are 
optimal with respect to erasures of the reconstruction system coefficients involved in 
the decoding process. We consider two different measures of the reconstruction error in 
a blind reconstruction algorithm. We also study the projective reconstruction system 
that best approximate an arbitrary reconstruction system, based on some well known 
results in matrix theory. Finally, we present a family of examples in which the prob- 
lem of existence of a dual projective system of a reconstruction system of this type is 
considered. 
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1 Introduction 



Many researchers have recently studied the problem of designing finite frames for the re- 
construction of a signal x in the absence of a number of (missing or delayed) coefficients 
(see |H [5j El [7J HU [121 E3 US])- m this context, the goal is to construct finite frames which 
minimize the (blind) reconstruction error of a signal when a number k of frame coefficients 
are "erased". 

In [2], the authors study the following problem: given a fixed frame T = {fi}^ for 
C d , find the alternate dual frames that minimize a measure of the (blind) reconstruction 
error when k frame coefficients are erased. Thus, unlike the previous setting, the interest 
is centered in a subset of the set of alternate dual frames of a fixed frame. As they notice, 
the canonical dual is not always the best choice despite the fact that it is associated to 
the Moore-Penrose pseudoinverse. Nevertheless, one of the main results in [14] exhibits 
conditions on the fixed frame J 7 which assure that the only optimal dual of T - in the sense 
previously described - is its canonical dual. 

Our research is in the vein of [H] but our interest is to set it in the context of fusion frames. 
However, since the definition of fusion frames depends strongly on the subspaces of the frame, 
the concept of dual fusion frame is not easy to handle. Therefore, we identify fusion frames 
with a subset of a broader class, namely the Reconstruction Systems (RS) (finite dimensional 
g-frames, using the terminology of W. Sun [2D]). Hence we fix a projective RS V (i.e. the 
analogue of fusion frames in the context of RS's), and we search for conditions that allow 
us to describe the dual RS's of V that minimize some measure of the (blind) reconstruction 
error. 

As in [TJ], our results show the existence of a unique dual RS (which is not necessarily 
projective) that is optimal for the erasure of 1 packet of coefficients. As in previous works in 
the subject, the error is measured in terms of the Frobenius norms of the so-called "packet 
lost operators" (the equivalent notion of the error operators in [U [11] E]). We present two 
different ways to perform this measure: the worst case error (WCE), i.e. the largest packet 
lost operator norm among all possible erasures, and the 2-error which is the euclidean norm 
of the vector of all packet lost operator norms. Each of these errors bound the norm of the 
error when a vector x is reconstructed with a blind reconstruction strategy, assuming that a 
frame coefficient is missing. 

For the WCE, following the approach of Han and Lopez, we impose some conditions on 
the norms involving the canonical dual elements to assure that the canonical dual is the 
only optimal alternate dual for one erasure. For the 2-error we show that, for projective 
RS's, there is always a unique optimal dual. In both cases, our results include a detailed 
description of the optimal dual RS. 

The fact that the optimal alternate duals of the fixed projective RS V considered above 
need not be projective motivates the study of best projective approximations of an arbitrary 
RS. We consider this problem in terms of well known results in matrix theory. 

On the other hand, assuming that the erased packets of coefficients correspond to a known 
subset of indexes of the RS V, we find conditions which assure that the RS corresponding 
to the remaining set of indexes is also a RS. Notice that this last fact allows to have perfect 
reconstruction with suitable dual RS's. We have included a section of examples in which we 
also consider the problem of the existence of projective duals of a given projective RS. 

In order to put our results in perspective and to keep the text self contained we have 
reviewed some known results concerning the structure of optimal duals in the context of 
fusion frames and ^-frames, including short proofs in some cases. It is worth pointing out 
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that some results in g-frames theory can be obtained from their analogues in frame theory, 
using the fact that g-frames can be considered as groupings of vectors of a frame. We have 
used this kind of argument whenever it was possible. Nevertheless, when an extra structure 
is imposed to the RS's new difficulties may arise; such is the case of projective RS's (i.e. our 
analogue of fusion frames), where there are orthogonal relations between subsets of vectors, 
or injective RS's. Indeed, as far as we know Theorem [7] is new even in the vector frame 
setting. In addition, by the nature of its formulation, the approximation by projective RS's 
and the existence of projective dual RS's of a fixed projective RS are intrinsic problems to 
the RS's setting. 

The paper is organized as follows: In Section [2] we recall some basic facts about general 
reconstruction systems and we fix some of the terminology used throughout the paper. In 
Section [3] we study the optimal dual systems for erasures of vector coefficients. In Section H] 
we consider the problem of perfect reconstruction when packets corresponding to a fixed set 
of indices of a RS V are lost. In SectionOwe describe the projective RS which is nearest to a 
fixed RS. In Section E] we give several examples of group RS's, and we consider the problems 
of existence of projective dual RS's of a given projective RS. 

Notations. 

Given H = C d and K = C n , we denote by L(7i, JC) the space of linear operators T : H — > /C. 
Given an operator T G L("H, /C), R(T) C K, denotes the image of T, kerT C % the null 
space of T and T* G L(/C,"H) the adjoint of T. If d < n we say that U G L("H,/C) is an 
isometry if U*U = I-u . In this case, U* is called a coisometry. We denote by I(d , n) the set 
of all isometries in L(7/,/C). 

If /C = % we denote by L(J-L) = L(H , Ti), by Gl{%) the group of all invertible operators in 
L(H), by L(7i) + the cone of positive operators and by Gl(7i) + = Gl(7i) fl L(7i) + . 

If T G L(H), we denote by <r(T) the spectrum of T, by rk T the rank of T, and by trT the 
trace of T. Given m G N we denote by I m = {1, . . . , m} C N and 1 = l m G M. m denotes the 
vector with all its entries equal to 1. 

On the other hand, A^„ )m (C) denotes the space of complex n x m matrices. If n = m we 
write A^ n (C) = M. n n {C), Ql{n) the group of all invertible elements of Ai n (C), U(n) the 
group of unitary matrices, Ai n (C) + the set of positive semidefinite matrices, and Ql(n) + = 

M n (C) + ngi(n). 

If W C TL is a subspace we denote by Pw G L(7i) + the orthogonal projection onto W, i.e. 
R(Pw) = W and ker Pw = W^. For vectors on C n we shall use the euclidean norm, but for 
matrices T G A4 n (C), we shall use both 

1. The spectral norm IITII = llTlLr, = max llTxII. 

Il^ll=i 

2. The Frobenius norm ||T|| 2 = (tr T*Tfl 2 = ( J2 \ T ij? ) 1/2 - Thi s 

norm is induced by 

the inner product (A, B) = tr B* A , for A, B e M n (C). 
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2 Basic framework of reconstruction systems 

In what follows we consider (m, k, (i)-reconstruction systems (see for example [151 HH1 120])- 
which are more general linear systems than those considered in [3J HI El El QT] and [IE] , that 
also have an associated reconstruction algorithm. 

Definition 1 Let m, d G N and k = (ki , . . . , k m ) G N m . 

1. We denote by K = /C m> k = ©j e i m . Sometimes we shall write each direct 
summand by /Q = C kl . 

2. Given a space U = C d we denote by L(m,k,d) = ® i&m L(U, 1Q = L(U,K). A 
typical element of L(m, k, d) is a system V = {V}iei m such that each Vi G L("H , /Q). 

3. A family V = {Vi}i e i m G L(m, k, d) is an (to, k, (i)-reconstruction system (RS) for "H if 

i.e., if Sy is invertible and positive. Sy is called the RS operator of V. In this case, 
the m-tuple k = (ki , . . . , k m ) G N m satisfies that trk = ki > d. 

4. The bounds of V as a RS are the numbers Ay = A m ; n (S'v) and By = \\Sy\\ sp . Observe 
that 

Av II ^ II 2 < (Sy x ) x ) — J2iei m II ^ x ll 2 — II ^ II 2 for every x eH . 

5. We shall denote by ftS(m, k, d) the set of all (m, k, d)-RS's for U = C d . 

6. The system V = {Vi}j £ i m G 1ZS(m, k, cf) is said to be injective if V* G L(/Q , "H) is 
injective (equivalently, if Vi V* G Gl{Ki) ) for every z G I m . 

We shall denote by HSx(m, k, d) the set of all injective elements of lZS(m, k, d). 

7. The system V is said to be projective if there exists a sequence v = (fi)iei m G of 
positive numbers, the weights of V, such that 

Vi V* = v] I Ki , for every i G I m . 

In this case, the following properties hold: 

(a) The weights can be computed directly, since each vi = ||V^|| sp . 

(b) Each Vi = ViUi for a coisometry U t G L(H , Ki). Thus V*V = vf Pr{v*) G L{H) + 
for every i G I m . 

(c) Observe that in this case Sy = J2iei m v 1 Pr(v*) as m fusion frame theory. 
We shall denote by V(m, k, d) the set of all projective elements of lZS(m, k, d). 

8. The analysis operator of the system V is defined by 

T v : H^f K = i6lm Ki given by T v x = (Vi x , . . . , V m x) , for x <E H . 
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9. Its adjoint is called the synthesis operator of the system V, and it satisfies that 
TZ:K = jeIm Id^n is given by ( ( Vi ) iElm ) = £ ieIm V* Vi . 
Using the previous notations and definitions we have that Sy = Ty Ty . A 
Examples 2 (Vector and fusion frames as RS's) 

1. As it was mentioned earlier, RS's arise from usual vector frames by grouping together the 
elements of the frame. Therefore, it is natural to expect that in the case k = t m , the set 
lZS(m, k, d) can be identified with the set of m-vector frames for H = C d . 

Indeed, let T = {fi}iei m G H™. For i G I m consider Vi : H — > C given by Vi(x) = (x, fi) for 
every x G H. Let Vjr = {^}«ei m and notice that 

Thus J 7 is a frame for % if and only if Vj- G lZS(m, 1, d). Actually, TZS(m, 1, d) = V(m, 1, d) 
because every functional is a multiple of a coisometry. Moreover, T^ T : "H — > ©i e i m C = C m is 
the usual analysis operator of J 7 . On the other hand, it is clear that elements in lZS(m, 1, d) 
correspond to vector frames for %. 

2. Let M w = (wi , Ai)iei m be a fusion frame for H = C d , with weights Wi > and subspaces 
Mi ^% with dimA/i = fcj for every i G I m . Its fusion frame operator is 

i£lm 

(see [HO [12] for detailed expositions of fusion frames). Let Ui G L{7i,C kl ) be a coisometry 
such that ?7*i7j = Pj^. , for every i G I m . Therefore, the system V = {^}iei m = {tf; K}iei m 
satisfies that Sy = SV m G Gl('H) + . Hence V G P(m,k, d) is a projective RS associated to 
M w . Observe that V has the same weights as M w and it also satisfies that each Mi = R(V*). 

Conversely, given V = {^}iei m £ V(m, k, d), the sequence AC = ( ||^||, R(V*) ) i f is a 
fusion frame such that Sy = Sjt/ W • Nevertheless the correspondence is not one to one, since 
any system of coisometries {t/j}tei m with (ker Ui) 1 - = Mi produces the same fusion frame 
M w = (wi , Mi)i£\ m . This phenomenon is similar to the correspondence of vector frames with 
one dimensional fusion frames. A 

Remark 3 In what follows we list some properties and notations about RS's: 
1. Given V = {V} ielm G KS(m, k, d) with S v = £ ieIm V*V , then 

J2Sy- 1 V*V i = I n , and ^ V*Vi S v ~ x = I H . (1) 

Therefore, we obtain the reconstruction formulas 

x = 5 V 1 V* (V x) = v * v ^ s v 1 x ) for ever y x eH . 
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2. For every V = {Vi\i e \ m G lZS(m, k, d), we define the system 



called the canonical dual RS associated to V. By Eq. (CQ), we see that 

T^T v =Y J S v - 1 V*V i = I H and S v# = ]T S v 1 V* V, S v 1 = S v 1 . 

Next we generalize the notion of dual RS's: A 

Definition 4 Let V = {Vi} ieIm and W = {Wi} ie i m G KS(m, k, d). 

1. We say that W is a dual RS for V if T v = I n , or equivalently if 

x = W* Vi x for every x G % . 

2. We denote by D(V) = {W G ftS(m, k, d) :T^T V = I H }, the set of all dual RS's for 
a fixed V G KS{m, k, d). Observe that D{V) ± since V # e D(V). 

A 

Remark 5 Let V G TZS(m, k, <i). Then W G -D(V) if and only if its synthesis operator 
is a pseudo- inverse of Ty . Indeed, W G -D(V) T-^Ty = . Observe that the map 

7ZS(m,k,d) 3 W i— > is one to one. Thus, in the context of RS's each (m, k, <f)-RS has 
many duals that are (m, k, <i)-RS's. This is one of the advantages of the RS's setting. 

Moreover, the synthesis operator of the canonical dual V* corresponds to the Moore- 
Penrose pseudo- inverse of Ty . Indeed, notice that Ty = Ty Sy l Ty G L(/C) + , so that it 
is an orthogonal projection. From this point of view, V* has some optimal properties that 
come from the theory of pseudo-inverses. A 



3 Optimal erasures and errors 

In coding theory, and using our terminology, a signal x G C d is transmitted encoded via 
the RS V = {Vi}i e i m in the form of m packets {Vix}j g i m . Then, the loss (or delay) of a 
number of packets in transmission is treated as if the corresponding components in the RS 
were "erased". Assuming that a number of erasures occurred and taking into account the 
redundancy of the RS, the reconstruction of the signal could be attempted using a dual 
RS of V. The accuracy of this process can be quantified by the norm of an error operator. 
For example, it is shown in [TT] that uniform Parseval frames (together with their canonical 
duals) are optimal for 1 erasures in the sense that they minimize a norm of the error operator. 
The authors also prove that equiangular uniform Parseval frames are optimal for 2 erasures 
(in case that such frames exist). In the context of Parseval fusion frames, Bodmann ([!]) 
shows that uniform Parseval fusion frames (a protocol associated to a uniformly weighted 
projective resolution of the identity, using his terminology) are optimal for 1 erasure, when 
the dimensions of the subspaces are equal. In addition, he shows that the uniform Parseval 
protocols that are optimal for multiple erasures have subspaces satisfying the so-called equi- 
isoclinic condition. 
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Similar problems have also been considered in different contexts. For example, in [12] 
the authors describe the structure of the optimal frames for one or multiple erasures, but 
in a non-deterministic setting. Indeed, in this case a random vector x is estimated from its 
fusion frame measurements using the Linear Minimum Mean Square Error. 

As it was mentioned in the introduction, our aim is to set the problem studied in [H] in 
the context of reconstruction systems. That is, given a fixed projective RS V, the goal is to 
characterize optimal alternate duals for r erasures. As in [13] , our results imply the existence 
of a unique dual RS (which is not necessarily projective) that is optimal for the erasure of 
1 packet of coefficients. Notice that in this context, optimal alternate duals for r erasures 
among those duals which are optimal for r — 1 erasures (see [31 El El El El US]). Therefore, 
by the nature of our results, the optimal duals for r erasures coincide with optimal duals for 
1 erasure. 

In order to describe the different measures of the reconstruction error when an arbitrary 
packet of coefficients of the fixed RS is erased, we consider the following notions. Given 
j elm, let 

Mj e L{K) given by M 3 ( (^ e iJ = ( \ 3 (i) 'Vi) ieIm , 
where tj : I m — » {0, 1} denotes the characteristic function of the set {j} C I m . Similarly, 
we consider the packet-lost operator Lj = Mf m \^y = Ijc — Mj . 

Given V = {Vi}« e i m G 1ZS(m, k, d), we shall consider a "blind reconstruction" strategy in 
case that some coefficient is lost. That is, assuming that the encoded information T^x e JC 
(for some x G %) is altered according to the packet-lost operator Lj , our reconstructed 
vector will be x = Lj T v (x), where W = {Wi}j £ i m G D(V) is some dual RS for V. Then 
the reconstruction error will be x — x = x — Lj T v (x) = Mj T v (x) = W*Vj x. 

In this case, we will use the Frobenius norm || • || 2 to perform the measure of the operator 
W*Vj . It is worth to note that, unlike the vector frame case, W*Vj is not necessary a rank- 
one operator, so its operator (spectral) norm does not coincide with its Frobenius norm. We 
consider that this is a suitable norm to perform the generalization of the results in [T3] to the 
RS setting, besides the fact that Frobenius norm has nice geometrical properties. Consider 
the m-tuple 

^(V, W) = {\\I-T^LjT v \\ 2 ) j&m = {\\T^MjT v \\ 2 ) j&m = (||^m) i6lm e K m , 

Notice that we can bound uniformly the reconstruction error in terms of the entries of this 
vector for the erasure of 1 packet of coefficients (for all m possible choices). In what follows 
we shall consider two different measures of the reconstruction error based on Ei(V, W), 
namely the 2-error and the (normalized) worst-case error. 

3.1 Minimizing the 2-error 

Let V = {Vi}i £ i m G 1ZS(m, k, d) and let us denote by 

. (2) 

We are interested in the characterization of those W G D(V) such that ||i?i(V, W)|| 2 = 

(2) 

e\ (V). In other words, we define the set of 1-loss optimal dual RS's for V with respect to 
\\Ex{V, -)|| a as 

Df\v) = {W G D{V) : ||£i(V, W)\\ 2 = e! 2) (V)} . 



ej 2) (V) = inf WE^V, W)|L = inf 
vveD(V) weD(v) 
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One of the advantages of this measure of the error is that for a projective V, without any 
further assumption, there is always only one dual RS in D^\v), and that this system can 
be explicitely computed (see Theorem [7] and Remark [8]). In particular, using the hierarchies 
previously described, this dual will also be optimal for r erasures. 

Remark 6 Fix V = {^};en m G V(m, k, d) and W = {W;}iei m eJ)(V). We give a matricial 
interpretation of the 2-error ||£'i(V, W)|| 2 , for the sake of clarity. Recall that V, W G 
L(m, k, d) = © igIm L(H , /Q). Hence they can be viewed as block m x 1 column vectors. 
For example, this allows to obtain the equality / = Ty = J2iei m Vi ; unfortunately 
this identity does not allow us to compute the value of ||i?i(V, W)|| 2 . 

On the other hand, the product T v = . gI Vi W* G L(fC) is an oblique projector, and 

\\TvTM\l = Ei,j& m \\ViWf\\l = Ei, j&m W*V4l, since the blocks V,W* G L{K 3 , /Q) are 
the different "entries" of the m x m block matrix Ty . 

The norm ||TyT^y|| 2 is minimal for W = V* (since TyT^ # is selfadjoint). Nevertheless, 
the square of the 2-error ||£i(V, ~W)\\ 2 = 11^7 Vi|| 2 * s ^ ne sc L uare °f the norm of 

the pinching matrix (i.e., block diagonal truncation) of TyT^,. As we shall see below, its 
minimum among all W G D{V) is not always attained at the Moore Penrose pseudo inverse 



T* # (see Remark HUD. A 

Theorem 7 Let V = {V}iei m G V(m,k,d) with weights v = (vi) ie i m . Then Df\v) = 
{Wo} i.e., there is a unique 1-loss optimal dual RS Wo for the 2-error. Moreover, if 

1. D G T(/C) is the block diagonal matrix D = @ v^ 2 I/c i , and 

iei m 

2. S VyD = T*DT V = £ P R{yn G GZ(ft)+ (sznce S VjD > (min v~ 2 ) -S v >0), 

then the optimal system is Wo = {v^ 2 Vi S'y^j-iei™ ■ In particular, T Wo = DT V S v l D . 

Proof. First we check that W G D(V). Indeed, T^ Q T V = S V ^ D T$DT V = S^ d S v ,d = I- 
Denote by P>i = Vi , the i-th. entry of Wo , for every i G I m . Consider a dual system 
W = {^}, eIm G D(V). Since each ^V? = ^ , then 

\\w*Vi\\ 2 2 = ti(y?WiW;vi) = triw^vv*) = v 2 tiiw.w*) = v 2 \\w*\\ 2 2 . (3) 



In particular, we can compute the 2-error for Wo : 

ll^(v,Wo)ii 2 = $:i|5;^ii^ i ^ = «r 2 W D « 

iei m ieim *ei m 

On the other hand, for every i 6 I m we have 



vtwnl = v 2 \\b* + {w; - B*)\\l 



,-2,10-1 



1^^*11? + ^ \\W* -B*\\ 2 2 +2Re (^ 2 tr [(W? - fl?)^] ) . (4) 
Let U = tr [ (W? - S?)S<] = ^ 2 tr [ W - 5 D V i s v!d] • Then we have that 

^^t, = tr( ^(WT-B^ViS^) = tr[(T^-T^)T v S^ D ] =0. 



ieim iei 
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since both W and V* are dual RS's for V. Therefore, summing over I m , the third summand 
of (@J vanishes and 

\\Ei(v, w)n 2 = e wmii i e % 2 iwii; 

1 £ ^II^d^II^ + ^IIW-^II 2 (5) 

*6l m 

.211 d*||2 © II 77 A) l/W MI2 



> E ^||£*|| 2 = ||£i(v, w ) 



Therefore W G ^ 2) (V). Moreover, if we take another W = {Wi} ieIm G D^(V), then Eq. 
d5J implies that — -B*|| 2 = for every z G I m , so that W = W . □ 

We say that a system V = {^}«ei m G k, d) is an uniform projective RS if the weights 
v = (vi)i e i m of V satisfy that v = v 1 for some t> > 0. 

Remark 8 The unique 1-loss optimal dual RS of Theorem [3, denoted by Wo G D\ (V), 
can be described in the following way: If v = (fj)j e i m are the weights of the system V G 
V(m, k, d), consider W = {£/i}iei m = {^}iei m > which is a uniform projective RS. Then 

Sv,d=J2 P m) = Su , V = v ■ U = {Vi Ui} i&m and W = v" 1 • U* , (6) 

ie Im 

because W = fa" 2 ^ ^v.^lieim = i v i l u i S^heim ■ A 

Corollary 9 Lei V = {Vijieim G V(m, k, d). Assume that V is uniform. Then the unique 
1-loss optimal dual RS for the 2-error is the canonical dual V*. 

Proof. If the weights of V are v = v 1 then, with the notations of Remark [HI we have that 
V = vU, S v = v 2 S u and V # = v~ l U* . Then we apply Eq. ® . □ 

Remark 10 For most systems V = {^}iei m G V(m, k, d) which are not uniform, the unique 
1-loss optimal dual RS for the 2-error Wo G D(V) obtained in Theorem [7] does not coincide 
with the canonical dual V*. For example, it is easy to see that Wo 7^ V* if all the weights are 
different and if the subspaces R(V*) are not mutually orthogonal. In particular, Wo 7^ V* 
whenever the weights are different and V has non zero redundancy. A 

3.2 Minimizing the worst-case reconstruction error 

Let V = {Vi}i e i m G V(m, k, d). For each W = {Wi}j £ i m G D(V) , we introduce the worst- 
case reconstruction error when one packet is lost with respect to the Frobenius norm. Now, 
we measure the error vector Ei(V, W) with the maximum of its entries, instead of the 
euclidean norm used in the previous section. 



ei (V) = inf \mV,W)\L = inf m f x\\T^M t T v \\ 2 = inf max \\W* . (7) 



9 



We define the set of 1-loss optimal dual RS's for V with respect to H-E^V, -JH^ as 

D 1 (V) ^ {WED(V) : WEkQt, w)|L = ei (v)} . 

The study of -Di(V) has been considered by Han and Lopez in [H] in the particular case 
of (m, l,d)-RS's for "H, i.e. usual vector frames. Indeed, since in such case the operators 
i G I m , T^MjTy are rank one operators, their Frobenius and spectral norms coincide, so the 
WCE coincide with the measure of the error used in [Hj. The use of the Frobenius norm 
in the definition of the WCE allows to extend naturally the results in [TJ] to the RS setting 
using similar techniques. 

Remark 11 Recall that KS(m, k, d) C L(m, k, d) = ieIm L(H , JQ = L(H, K). If we fix 
an injective V = {Vi}i e i m G IZSxijn, k, d), then the map || • ||v : L(m, k, d) — > R + given by 

||W|| V = \\Ei(V, W)L = max||W^|| 2 for W = {Wi} i6lm G L(m, k, d) 

is a norm in L(m, k, d). Indeed, the only non trivial condition is the faithfulness. But the 
fact that V G IZSxijn, k, d) (i.e. is surjective for every i G I m ) assures that ||W||v = ==>■ 
\\W* Vi\\ 2 = for every % G l m W = 0. 

Since D(V) is closed in L(m,k, d) with the usual norm and all norms are equivalent on 
L(m, k, d), then D(V) is || ■ ||y -closed in L(m, k, d). Therefore there exist elements W G -D(V) 
such that 1 1 W || v — mm ll-^llv = e i("^)- Indeed, the intersection of D(V) with a fixed 

closed ball is a compact set. On the other hand, D(V) is convex (actually it is an affine 
manifold). Since every norm is a convex map, we have proved the following result: A 

Proposition 12 Let V G lZSx(m, k, d) be an injective system. Then the set -Di(V) of 
1-loss optimal dual RS's for V is non-empty, compact and convex. □ 

Theorem 13 Let V = {^}iei m G "P(m, k, d) weights v = (fi)i 6 i m . 1/ 

ll-Sv^Vr^lla = "ill^v^WjIla = c /or every z G I m , 

t/ien V* ; £/ie canonical dual RS ofV, is the unique 1-loss optimal dual RS for V (and hence 
the r-loss optimal dual RS for every r). In other words, -Di(V) = {V*}. 

Proof. By Proposition [T2"| there exists some W = {W^}j e i m G -Di(V). Then 

|| W||v =max||^\/ i || 2 < rnaxHSv 1 ^* V5|| a = l|V # ||v = c . 

If we denote each ViS^ 1 = C i} then HW/^I^ < c = ||C*V^||^ for every % G I m . Recall 
that, by Eq. ©, 1 1 W^*T^ 1 1 ^ = u?||W7|||, since V-V? = v 2 I K . . Similarly, we get that each 
\\C*Vi\\l = v? \\C*\\ 2 2 . Therefore \\Wi\\ 2 2 < \\C*\\ 2 2 for every % G I m . Note that 

\\W*\\l = \\C* + (W* -C*)\\ 2 2 = ||C*||2 + ||V^*-C*||2 + 2Re(tr [{W* - C*)d] 
and hence ||W^ - C*\\ 2 2 + 2 Re ( tr [ (W* - C*)Cj ) < 0, for every i G l m . Finally, 

tr [ (W* - C*)C t ] = tr [ (T^ - T^)T V S^ 
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since both W and V* are dual RS's for V. Then 

o < Yl w w * - c *wi = E iw - c *wi + E 2 Re ( tr - c ^ ) ^ - 

«ei m «ei m iei m 

which implies that W = {Wi} ieIm = {a}iei m = V # . □ 

A system V G 1ZS (m, k, <i) is called a protocol for 7/ if Sy — 1% . This notion appears in 
[I], PS] (see also j5], where protocols are related to C*-encodings with noiseless subsystems). 

Corollary 14 Let V = {^}iei m £ "P(m, k, d) fre a projective protocol for 7-L (i.e. Sy — I) 
such that \\V*Vi\\ 2 = vf k] /2 = c for every iel m . Then D X {V) = {V*} = {V}. 

Proof By hypothesis Sy = In , an d hence H^y 1 V^*Vi|| 2 = ||V^*T^|| 2 = c for every i e I m . 
Thus, the previous theorem can be applied in this case. 

Remark 15 Examples of projective protocols as in Corollary [TH are the equi- dimensional 
uniform projective protocols i.e., {^}«ei m £ V(m,kt,d) that are uniform. These are the 
analogues of the so-called uniform fusion frames. 

The (m, 1, d) case of Theorem IT51 is a rephrasing of [HI Thm 2.6], since all vector frames are 
projective as RS's. Moreover, we can conclude from the examples P31 Section 3] that the 
optimal dual system W G -Di(V) may not be the canonical dual RS and may be not unique 
for a general V G V(m, h,d). A 



4 Stability of RS's under erasures of coefficient packets 

In this section we consider a different approach to the erasures problem. Indeed given a 
fixed RS V, assume that we can identify the set J C I m such that {Vi are the missing 
or delayed packets of coefficients in the transmission of the signal x G H. In this case 
we shall state conditions which assure that the system Vj = {Vi\i & i m \j corresponding to 
the remaining set of indexes is still a RS. Notice that this last fact allows to have perfect 
reconstruction with suitable dual RS's. Hence, an explicit computation of the canonical dual 
(Vj)* is given in this case. 

The following statement is a generalization and a slight improvement of similar results of P. 
Casazza and G. Kutyniok [7] and Asgari pQ for fusion frames. 

Proposition 16 Let V = {^}tei m G 7ZS(m, k, d) with bounds Ay , By. Fix a subset J C I m 
and consider the matrix Mj = Id — J2i<=j Sy 1 e A^ n (C). Then, 

Vj = (VA mm \j is a RS for H = C d <=>• Mj G Ql(d) . (8) 

In this case we can compute the following data for Vj : 

1. The frame operator Syj = Mj Sy . 

2. The bounds of Vj can be estimated by ^ < Ay, and By, < By . 

\\Mj || 

3. The canonical dual can be characterized as 

(Vj)* = {Vi S^Uj = {Vf Mj% i3 d d ( V #) j ■ MJ 1 . 
That is, (Vj)* is the truncation of the canonical dual V* modified with Mj 1 . 
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Proof. It is straightforward to check that Mj = Syj S v =>- Syj = Mj Sy . This last fact 
implies the equivalence of Eq. (151) . On the other hand, 

Ay \\MJ X \\- 1 < UMjSy)- 1 ]]- 1 = \\SyX 1 = A v.j ■ 

The fact that < Syj < Sy assures that By } < By. Let us denote V # = {W^}^^ and 
= {Wi}i^j . Then the formula Sy a = Mj Sy gives the equality 

(V#)j ■ MJ 1 d = f {W l Mj l } HJ = {Vi Sy 1 Mj l } HJ = {Vi S^htj = (Vj)# . □ 

Remark 17 In [7j P. Casazza and G. Kutyniok stated the sufficient (and easily computable) 
condition ^ ||Vi||^ p < Ay for the invertibility of the matrix Mj of Proposition [T6l Indeed, 

if II S ^7 Vi\\sp < Ay (compare with the condition ^ ||^||L < Ay), then || I d — Mj \\ sp < 1 

i&J ieJ 

and Mj e Ql(d). They also give the estimation Ay — \\Vi \\% < Ay T . This follows from 

PropfTBI 

A-v-YnzAVitp < Av-\\YKjViVi\U< pit < a Vj . 

The result of Asgari in pQ is similar to Eq. flSJ), but stated for fusion frames and assuming 
that | J | = 1, with a different lower bound for Vj . A 



5 Approximation by projective RS's 

Given a fixed V G V(m, k, rf), notice that the optimal dual RS's obtained in Theorems [TJ and 
[131 , are not projective RS's, in general. Although there could be some projective elements 
in D(V) (we shall focus this problem in the following section), we are interested in those 
(m, k, (i)-projective RS's that are closest, with respect to some distance function, to a fixed 
S = {Si}i e i m G D(V) which has some desired properties. Given W G lZS(m, k, d), we 
consider 

d(S,W) = \\T S -T W \\ 2 = \\T*-T^\\ 2 , 

the distance between their synthesis (or analysis) operators. Hence, we seek for Wo G 
V(m, k, d) that minimize d(S , W) among the projective RS's. In what follows we will 
describe the structure of such (unique) minimizers in case S is an injective RS. As one would 
expect, its "directions" are the coisometries of the polar decompositions of the coordinate 
operators Si of S, while its weights are the "averages" of their singular values. We need first 
some preliminary results: 

Given k, n G N such that k < n, we denote by 

X(k ,n) = {Ue L{C k , C n ) : U*U = I k } , 
the set of isometries from C k into C n . The following result can be found in [2]: 

Lemma 18 Let k,n G N such that k < n, and let A G A4 ni fc(C) be a full rank matrix with 
polar decomposition A — U \A\ with U G X{k , n). Then \\A — U\\ 2 = min \\A — V\\ 2 . 

Vel(k,n) 

□ 
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Recall that for every A G Ai n> k(C) its polar decomposition satisfies 

A = U\A\ = \A*\U U*A = \A\ and A* = U*\A*\ , (9) 

where U G .M n> fc(C) has kerU = ker A. As expected, the aproximation of a RS by a projec- 
tive RS {<Si}iei m relies on the coisometries that best approximate each £», which are deter- 
mined by the polar decomposition of S{. Nevertheless, one has to determine the appropiate 
weights which combine those coisometries. That is done in the following Proposition. 

Proposition 19 Let S = £ TZSz(m, k, d). Then there exists a unique 

W G V(m,k,d) such that d(S , W ) = min d(W , S) , (10) 

tr\S\ 

and it is given by Wo = £^}iei m where each cti = — - — — and Si = Ui\Si\ is the polar 
decomposition of each Si . 

Proof. Let W G V(m, k, d) be a system such that the minimum in (fIU|) is attained at W. 
Denote by w = (wi) ie i m G R+ the weights of W. Notice that 

m 

l|35-2wli; = 53l|Sr-W7li; and each WW = w? I* , 

Thus, each isometry w^ 1 W* G X{ki , d) attains the minimum in the optimization problem 

c* c* 

\\w- l Wt- — \\ = min ||X-^|| , 

11 Wi " 2 xei{ki,d) 11 " 2 

where, by hypothesis, each w~ l S* is a full rank linear transformation. By Lemma [TS] we get 
that w~ l W* = U* , the isometry of the polar decomposition 

w^S* = U* {w; 1 S* | = wi 1 ( U* \S* | ) M> Si = Ui\Si\ and ker = ker Si , 
and hence Wi — WiUi . Next we show that each Wi = tr [ s ^ . Fix i G I m . Then 

- Wit/j|| 2 = min \\Si - a UAL . 

Therefore Wi ■ \\Ui\\ 2 is the norm of the orthogonal projection of Si to the line Kf/j , using 
the IR-inner product (A, B) = Re [tr(B*A)]. It can be computed explicitly: 



< tr \Sj\ j| tr 

1 1 1 1 2 1 1 1 1 2 



s K 



u « Ha 



PW3)ll» = «>i -11^ 



* II 2 5 



tr I Si I 

for every z G I m . Then we obtain the equalities Wi = \\Ui\\~ 2 tr\Si\ = — — — . □ 
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6 Examples 



In this section we present a variety of examples related with the previous sections. First, 
we exhibit a family of (projective) RS's which satisfies the hypotheses of Corollary |9] and 
Theorem [TBI the so-called group RS. In particular, the canonical dual of a group RS (which 
is also a group RS) is the unique 1-loss optimal dual RS for the 2-error as well as for the 
WCRE. The remaining subsections are devoted to the study of particular cases of RS where 
D(V) n V(m, k, d) ^ . The first examples show that for certain projective systems we can 
explicitly construct projective duals, which in general will not coincide with the canonical 
duals. The last example describes a Riesz RS whose unique dual RS (i.e. the canonical dual) 
is not projective. This leads to a characterization for Riesz RS's with projective canonical 
dual. 

6.1 Group reconstruction systems 

We begin by rephrasing some basic notions and results from [13] in the RS's setting: Let 
/Co = C fc and V £ L(7-L , /Co)- Given a unitary representation G 3 g >->■ U g £ U{%) of a finite 
group G in U (H) we say that 

V(G,V) d = f {VU g } geG 

is a G-reconstruction system (shortly, G-RS) if V(G , V) £ 1ZS(m, kt, d), where m = \G\. 
In this case the space K = /C™ . If S = V*(JC ), this is equivalent to the fact that 

spanj |J U g {S)}=H, 

where V £ L(H , /C) is the base operator for V(G , V). This definition of G-RS reduces 
to that of G-frame in the vector frames setting. Following [13] and [10] we state a series of 
properties of G-RS's whose proofs are similar to the frame case: Fix V £ L("H , /Co). 

1. Observe that the system V(G , V) is: 

(a) Projective (and uniform) if in addition VV* = v 2 I/c for some v > ; 

(b) Injective if in addition V* is injective, in which case also R(V U g ) = K, for every 
geG. 

2. Notice that the RS-operator of V(G , V) has the following structure: 

s G ,v = s v[G , v) = Y J u;v*vu g = Y J Ug-^v*vu g . 

geG g&G 

3. The RS-operator Sq, v (and therefore S^y) commutes with the unitary representation 
of G: 

Uh ■ Sg,v = S GtV ■ U h for every heG. (11) 

4. In particular, the canonical dual of a G-RS is another G-RS: 

V(G , V)* = {V U g Sa] v } 9 e a = {V S G ] V U g } ge G = V(G,V S G ] V ) . 
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In order to apply our previous results, assume now that the base operator V G L(7i, /C) 
satisfies VV* = v 2 J^. Then V(G, V) is an equi- dimensional uniform projective RS. Hence, 
Corollary [9] implies that the 1-loss optimal dual RS for V(G , V) for the 2-error is its canonical 
dual V(G , V)*. In this case then V(G , V) also satisfies the hypothesis of Theorem [T3j 
V(G , V) is a G-projective RS, and 

U0s£ v t/; tg|| 2 = || u; (s G \ v*v) u g \\ 2 = \\s G ] v v*v\\ 2 = c . 

Thus, the canonical dual of such G-RS's is the unique 1-loss optimal dual for the worst-case 
error. If the base operator V G LiTi, /C) is surjective, so that V(G , V) is an injective RS 
then, using Proposition [191 the projective RS nearest to V(G , V)* can be computed in the 
following way: For every g G G, we have that 

\vs G ] v u g \ 2 = u; s G \ v* v stfy u g = u; \vs G ] v \ 2 u g . 

Taking square roots at both sides, we get that \VU g S G v \ = U*\VS G v \U g for every g G G. 
Therefore, if we consider the polar decomposition V S^y = W\V S G v \ of V S^y , then 
also 

vs G ) v u g = {wu g ) {u;\vs G ] v \u g ) = {wu g ) \vu g s G ] v \ 

is the polar decomposition of each entry V S^y U g of V(G , V)*. In conclusion, if we denote 

tr IVS^VI 

w = °' vl , then V(G , w W) = { w W U g } geG 

is the best projective approximation of V(G , V)* provided by Proposition [19j It is clear 
from the previous computations that it is again a G-RS. 



6.2 Dual projective systems 

Next we consider an example of a system V = {Vi} ie i m G V(m, k, d) with projective dual 
systems but such that V # ^ V(m,k,d) : 

Example 20 Let d = 3, m = 2 and k = (2, 2). Let V x and V 2 G L(C 3 , C 2 ) be given by 

Vi(x, y, z) = (y, z) and V 2 (x,y, z) = (x, z) for every (x, y,z)EC 3 . 
Then V = (Vi , V£) G P(m, k, d) with weights 1 2 . If S x = {e^ and <S 2 = {62}^, then 

/ 1 

S v = V* V x + V* V 2 = P Sl + P S2 =[ 1 

\ 2 

Therefore V # ^ V(m,k,d) since 5y 1 V^*(m, v ) = (0,-u, |) for (u,v) G C 2 , so that the entry 
Vi S^ 1 of V* is not a multiple of a coisometry. 

Let W = (Wi , W 2 ) G P(m, k, d) and assume that T^T V = W x * V x + V 2 = h ■ Denote by 
{v\ , v 2 } the canonical basis of C 2 and by {ei , e 2 , 63} that of C 3 . Then, easy computations 
using the definition of V show that 

e 3 = W* v 2 + W; v 2 , e 2 = W*V t e 2 = W* v± and e x = W*V 2 e x = W* v x . (12) 




15 



The last two equalities show that both W* and W£ should be isometries with weight 1. But 

2 ^ 1 2 



1 . /o 

in this case ||W* i^H = 11^2 — 1 an d their sum also has norm one. Let oo = t> + i^r ■ 



Then |w| = 1 but u + 57 = 1. Then we can define , W 2 * G L(C 2 , C 3 ) by 

Wi(x,y) = (0,x,uy) and W^ix.y) = (x, 0,57?/) for every (x,y) G C 2 . 

The W* are isometries and satisfy the three conditions of (TT2|) . Therefore, the system 
W = (Wi , W 2 ) lies in P(m, k, d) and it is a dual-RS for V. 

Nevertheless, if we consider V\ and V2 as operators in L(IR 3 , R 2 ), then such a W can not 
exists in the real case. Indeed, looking at Eq. (fl2l) . we can deduce that W* V2 G ^j -1 and 
W 2 * ^2 G {ei}- 1 . These facts, together with the equality e% = W* v% + W 2 * t>2 imply that both 
W* v<2 and W 2 * v 2 G span{e 3 }, which is impossible in the real case. A 

Example 21 We can generalize Example [2D] in the following way: Assume that the projec- 
tive system V = {Vi} ie i m G V(m, k, d) has the property that all the projections p = Pr(v*) 
are pairwise commuting. In this case D(V) fl V(m, k, d) . 

Indeed, suppose first that all the weights of V are 1. Then Sy = X^en ^ • The commutation 
hypothesis assures that, by taking all the possible intersections among the ranges of the 
projections p , we get a family of projections (Qj)jei n such that 

1. Qi Qj = if i ^ j and J2j E i n Qj = l u ■ 

2. S v = Y^jeK r i with r i e lm for ever y 3 e ln ■ 

3. For every i G I m there exists Jj C I n such that P = ^- 6J . • 

We construct the system W = {Wi} iElm G P(V) n V(m,k,d) as follows: let Wi = Villi, 
where 

1 y/3 

Ui = £y Qj for some G { 1 , — 1 , u , 57 } , where u — - + i —— . 

Note that since all = 1, then W G P(m, k, d). A careful selection of these coefficients, 
taking account the parity of the numbers rj = where Sj = {i G I m : j G Jj} for 

j G I„ , allows to find such a W such that W G -D(V). The general case follows from 
the previous case. Indeed, if V has weights v = (fj)j e i m , we replace the previous W by 

w v = K 2 w*W " a 

Remark 22 (Projective dual pairs) The previous example gives a method to construct 
pairs of projective RS's (V, W) such that W G -D(V). Moreover, this method shows that for 
every choice of parameters (m, k, d) such that ^2 ie j m K > d there exist V, W G V(m, k, d) 
such that W G -D(V). Indeed, to find such a projective dual pair (V, W) we construct 
V G P(m, k, d) in such a way that the projections p = V^*T^ for i 6 I m are pairwise 
commuting (i.e. that are simultaneously diagonalizable by an orthonormal basis of %). 
Then, we can apply the construction in the example above to obtain explicitly the projective 
dual W. 

These facts show that projective dual pairs are indeed more frequent than projective pro- 
tocols i.e. RS's V such that (V, V) is a projective dual pair, since it is known that there 
are choices of parameters (m, k, d) for which no projective (m, k, d)- protocols exist (see [XT] 
Example 3.1.2.]). A 
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6.3 Riesz reconstruction systems 

The elements of lZS(m, k, d) are called Riesz RS's if dim/C = trk = d = dimH . In this 
case, every V = {Vi}i e i m G 1ZS(m, k, d) has the property that both the synthesis operator 
Ty and the analysis operator Ty are invertible. Also T* # = 1 Ty = T^ 1 , and 

T-L = R(V*) © • • ■ © R{Vm) (direct sum, but not necessarily orthogonal) , (13) 

because the sum gives always % but in this case the sum must be direct by dimensional 
reasons. On the other hand, if V = {^}j S n m G TZS(m, k, d) is a Riesz RS for 7-L then 

D(V) = {V*} , (14) 

since the only left inverse of T v is T^ 1 = T^ # . Recall that T v = 7« for every W G -D(V) , 
and that the map lZS(m, k, d) 9 W (->■ T-^ is one to one. 

Example 23 Let d = 4, m = 2 and k = (2,2). We now construct a (necessarily) Riesz 
V G P(m, k, d) such that D(V) H P(m, k, d) = 0. Let 14 , V 2 G 7(C 4 , C 2 ) be given by 

Vi(xi , x 2 , x 3 , x 4 ) =(a;i,x 2 ) and ^(^i , x 2 , x 3 , x 4 ) = {x 3 , £ ^ i ) , (15) 

for (xi , x 2 , x 3 , x±) G C 4 . It is easy to see that V = (Vi , V 2 ) G V(m, k, d) with weights 
(1, 1). Let us denote by S = ker V 2 = spanjei , e 2 + e 4 } C C 4 . Given W G -D(V), the equality 

W* V x + Vr 2 * V 2 = I H 

implies that W{ V\ x — x for every x G S. Then W* G L(C 2 , C 4 ) is completely determined 
as the inverse of V% L : S 1 — )■ C 2 . But we have that 

||y iei || = |NI=l while ||li(e 2 + e 4 )|| = ||e 2 || = ^ ^ ■ 
Then V\\ s is not a multiple of an isometry and neither is . 

We can enlarge the previous example in order to get a RS with redundancy and without 
projective duals. Indeed, consider Vo = (Vi , V 2 , V 3 ) G V(3 , (2, 2, 2) , 4) obtained from V by 
adding any coisometry V3 G L(C 4 , C 2 ) such that also ker V3 = S. Then, arguing as before, 
we conclude that there is no W = (Wi , W 2 , W 3 ) G D(Vq) such that PVi is a multiple of a 
coisometry. A 

Remark 24 Assume that V = {V^}i e i m G lZS(m, k, d) is a Riesz RS. Then, arguing as in 
the previous example, it is easy to see that the following conditions are equivalent: 

1. D(V) nV(m,k,d) ^ 0. 

2. V* G V(m,k,d). 



3. If we denote by Si — f] ker Vi = ( R(V*) ) , then V^| s G 7(5*4 , ^Q) is a multiple of 
an isometry, for every i G I m . 

These conditions are fulfilled if the sum of Eq. fflBl is orthogonal. Also if every ki = 1. But 
there exist other cases in which D(V)nV(m, k, d) 7^ 0. For example, if we take the operator 
Vl of Eq. ( TlBl) , and consider 

V 3 (xi , x 2 , x 3 , x 4 ) = {xi - x 3 , x 2 - x 4 ) for (xi , x 2 , x 3 , x 4 ) G C 4 , 

then the condition 3 is satisfied by V = (Vi , V 3 ) G V(2, (2, 2), 4). A 
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